Development of Clinical Decision Support Systems based on Mathematical Models of Physiological Systems

In the last years of research, I focused my studies on different physiological problems. Together with my supervisors, I developed/improved different mathematical models in order to create valid tools useful for a better understanding of important clinical issues. The aim of all this work is to develop tools for learning and understanding cardiac and cerebrovascular physiology as well as pathology, generating research questions and developing clinical decision support systems useful for intensive care unit patients. I. ICP-model Designed for Medical Education We developed a comprehensive cerebral blood flow and intracranial pressure model to simulate and study the complex interactions in cerebrovascular dynamics caused by multiple simultaneous alterations, including normal and abnormal functional states of auto-regulation of the brain. Individual published equations (derived from prior animal and human studies) were implemented into a comprehensive simulation program. Included in the normal physiological modelling was: intracranial pressure, cerebral blood flow, blood pressure, and carbon dioxide (CO2) partial pressure. We also added external and pathological perturbations, such as head up position and intracranial haemorrhage. The model performed clinically realistically given inputs of published traumatized patients, and cases encountered by clinicians. The pulsatile nature of the output graphics was easy for clinicians to interpret. The manoeuvres simulated include changes of basic physiological inputs (e.g. blood pressure, central venous pressure, CO2 tension, head up position, and respiratory effects on vascular pressures) as well as pathological inputs (e.g. acute intracranial bleeding, and obstruction of cerebrospinal outflow). Based on the results, we believe the model would be useful to teach complex relationships of brain haemodynamics and study clinical research questions such as the optimal head-up position, the effects of intracranial haemorrhage on cerebral haemodynamics, as well as the best CO2 concentration to reach the optimal compromise between intracranial pressure and perfusion. We believe this model would be useful for both beginners and advanced learners. It could be used by practicing clinicians to model individual patients (entering the effects of needed clinical manipulations, and then running the model to test for optimal combinations of therapeutic manoeuvres). II. A Heterogeneous Cerebrovascular Mathematical Model Cerebrovascular pathologies are extremely complex, due to the multitude of factors acting simultaneously on cerebral haemodynamics. In this work, the mathematical model of cerebral haemodynamics and intracranial pressure dynamics, described in the point I, is extended to account for heterogeneity in cerebral blood flow. The model includes the Circle of Willis, six regional districts independently regulated by autoregulation and CO2 reactivity, distal cortical anastomoses, venous circulation, the cerebrospinal fluid circulation, and the intracranial pressure-volume relationship. Results agree with data in the literature and highlight the existence of a monotonic relationship between transient hyperemic response and the autoregulation gain. During unilateral internal carotid artery stenosis, local blood flow regulation is progressively lost in the ipsilateral territory with the presence of a steal phenomenon, while the anterior communicating artery plays the major role to redistribute the available blood flow. Conversely, distal collateral circulation plays a major role during unilateral occlusion of the middle cerebral artery. In conclusion, the model is able to reproduce several different pathological conditions characterized by heterogeneity in cerebrovascular haemodynamics and can not only explain generalized results in terms of physiological mechanisms involved, but also, by individualizing parameters, may represent a valuable tool to help with difficult clinical decisions. III. Effect of Cushing Response on Systemic Arterial Pressure. During cerebral hypoxic conditions, the sympathetic system causes an increase in arterial pressure (Cushing response), creating a link between the cerebral and the systemic circulation. This work investigates the complex relationships among cerebrovascular dynamics, intracranial pressure, Cushing response, and short-term systemic regulation, during plateau waves, by means of an original mathematical model. The model incorporates the pulsating heart, the pulmonary circulation and the systemic circulation, with an accurate description of the cerebral circulation and the intracranial pressure dynamics (same model as in the first paragraph). Various regulatory mechanisms are included: cerebral autoregulation, local blood flow control by oxygen (O2) and/or CO2 changes, sympathetic and vagal regulation of cardiovascular parameters by several reflex mechanisms (chemoreceptors, lung-stretch receptors, baroreceptors). The Cushing response has been described assuming a dramatic increase in sympathetic activity to vessels during a fall in brain O2 delivery. With this assumption, the model is able to simulate the cardiovascular effects experimentally observed when intracranial pressure is artificially elevated and maintained at constant level (arterial pressure increase and bradicardia). According to the model, these effects arise from the interaction between the Cushing response and the baroreflex response (secondary to arterial pressure increase). Then, patients with severe head injury have been simulated by reducing intracranial compliance and cerebrospinal fluid reabsorption. With these changes, oscillations with plateau waves developed. In these conditions, model results indicate that the Cushing response may have both positive effects, reducing the duration of the plateau phase via an increase in cerebral perfusion pressure, and negative effects, increasing the intracranial pressure plateau level, with a risk of greater compression of the cerebral vessels. This model may be of value to assist clinicians in finding the balance between clinical benefits of the Cushing response and its shortcomings. IV. Comprehensive Cardiopulmonary Simulation Model for the Analysis of Hypercapnic Respiratory Failure We developed a new comprehensive cardiopulmonary model that takes into account the mutual interactions between the cardiovascular and the respiratory systems along with their short-term regulatory mechanisms. The model includes the heart, systemic and pulmonary circulations, lung mechanics, gas exchange and transport equations, and cardio-ventilatory control. Results show good agreement with published patient data in case of normoxic and hyperoxic hypercapnia simulations. In particular, simulations predict a moderate increase in mean systemic arterial pressure and heart rate, with almost no change in cardiac output, paralleled by a relevant increase in minute ventilation, tidal volume and respiratory rate. The model can represent a valid tool for clinical practice and medical research, providing an alternative way to experience-based clinical decisions. In conclusion, models are not only capable of summarizing current knowledge, but also identifying missing knowledge. In the former case they can serve as training aids for teaching the operation of complex systems, especially if the model can be used to demonstrate the outcome of experiments. In the latter case they generate experiments to be performed to gather the missing data.

Ruolo delle alterazioni mitocondriali indotte dall'ipossia nella patogenesi della malattia di alzheimer

Le alterazioni della funzionalità mitocondriale detengono un ruolo cruciale nella patogenesi della malattia di Alzheimer (AD), sostenendo il processo neurodegenerativo attraverso meccanismi quali la riduzione della disponibilità energetica e la iperproduzione di ROS. Alle numerose ipotesi di patogenesi dell’AD, si è recentemente affiancata la cosiddetta ipotesi vascolare. Nei soggetti AD è stata riscontrata una significativa riduzione della disponibilità di ossigeno a livello neuronale (ipossia neuronale). Da numerosi studi è poi emerso che l’ipossia gioca un ruolo fondamentale nello sviluppo dell’AD contribuendo a più vie patogenetiche contemporaneamente. Tuttavia, non sono stati ancora chiariti tutti i meccanismi attraverso cui l’ipossia esplica la sua azione di danno. Lo scopo di questo studio è stato quello di contribuire a chiarire il ruolo patologico dell’ipossia nell’AD, analizzando principalmente le alterazioni della funzionalità mitocondriale indotte dalla riduzione della disponibilità di ossigeno. Nella prima fase dello studio cellule PC12 sono state coltivate in presenza di β-amiloide e ipossia. In questo modello abbiamo osservato un potenziamento dei fenomeni di deplezione dell’ATP e di generazione delle ROS indotti dalla Aβ quando anche l’ipossia era presente come fonte di danno cellulare, ipotizzando per i due fattori un effetto congiunto di tipo additivo. Nella seconda fase abbiamo esposto all’ipossia fibroblasti prelevati da pazienti AD portatori di mutazioni a carico dei geni APP e PSEN. La presenza di mutazioni predisponenti ad un fenotipo AD era in grado di determinare un danno bioenergetico e ossidativo. Le alterazioni bioenergetiche riscontrate in normossia risultavano ulteriormente potenziate quando i fibroblasti erano coltivati in ipossia, mentre lo stato di stress ossidativo veniva evidenziato solo in condizioni ipossiche. Sulla base dei risultati finora conseguiti si può ipotizzare che uno dei meccanismi attraverso cui l’ipossia esplica la sua azione di danno nella AD, possa essere dovuto alla capacità di potenziare ulteriormente le alterazioni della funzionalità mitocondriale.

Angiogenesis and angioregression gene expression analyses in swine corpus luteum

The corpus luteum (CL) lifespan is characterized by a rapid growth, differentiation and controlled regression of the luteal tissue, accompanied by an intense angiogenesis and angioregression. Indeed, the CL is one of the most highly vascularised tissue in the body with a proliferation rate of the endothelial cells 4- to 20-fold more intense than in some of the most malignant human tumours. This angiogenic process should be rigorously controlled to allow the repeated opportunities of fertilization. After a first period of rapid growth, the tissue becomes stably organized and prepares itself to switch to the phenotype required for its next apoptotic regression. In pregnant swine, the lifespan of the CLs must be extended to support embryonic and foetal development and vascularisation is necessary for the maintenance of luteal function. Among the molecules involved in the angiogenesis, Vascular Endothelial Growth Factor (VEGF) is the main regulator, promoting endothelial cells proliferation, differentiation and survival as well as vascular permeability and vessel lumen formation. During vascular invasion and apoptosis process, the remodelling of the extracellular matrix is essential for the correct evolution of the CL, particularly by the action of specific class of proteolytic enzymes known as matrix metalloproteinases (MMPs). Another important factor that plays a role in the processes of angiogenesis and angioregression during the CL formation and luteolysis is the isopeptide Endothelin-1 (ET-1), which is well-known to be a potent vasoconstrictor and mitogen for endothelial cells. The goal of the present thesis was to study the role and regulation of vascularisation in an adult vascular bed. For this purpose, using a precisely controlled in vivo model of swine CL development and regression, we determined the levels of expression of the members of VEGF system (VEGF total and specific isoforms; VEGF receptor-1, VEGFR-1; VEGF receptor-2, VEGFR-2) and ET- 1 system (ET-1; endothelin converting enzyme-1, ECE-1; endothelin receptor type A, ET-A) as well as the activity of the Ca++/Mg++-dependent endonucleases and gelatinases (MMP-2 and MMP-9). Three experiments were conducted to reach such objectives in CLs isolated from ovaries of cyclic, pregnant or fasted gilts. In the Experiment I, we evaluated the influence of acute fasting on VEGF production and VEGF, VEGFR-2, ET-1, ECE-1 and ET-A mRNA expressions in CLs collected on day 6 after ovulation (midluteal phase). The results indicated a down-regulation of VEGF, VEGFR-2, ET-1 and ECE-1 mRNA expression, although no change was observed for VEGF protein. Furthermore, we observed that fasting stimulated steroidogenesis by luteal cells. On the basis of the main effects of VEGF (stimulation of vessel growth and endothelial permeability) and ET-1 (stimulation of endothelial cell proliferation and vasoconstriction, as well as VEGF stimulation), we concluded that feed restriction possibly inhibited luteal vessel development. This could be, at least in part, compensated by a decrease of vasal tone due to a diminution of ET-1, thus ensuring an adequate blood flow and the production of steroids by the luteal cells. In the Experiment II, we investigated the relationship between VEGF, gelatinases and Ca++/Mg++-dependent endonucleases activities with the functional CL stage throughout the oestrous cycle and at pregnancy. The results demonstrated differential patterns of expression of those molecules in correspondence to the different phases of the oestrous cycle. Immediately after ovulation, VEGF mRNA/protein levels and MMP-9 activity are maximal. On days 5–14 after ovulation, VEGF expression and MMP-2 and -9 activities are at basal levels, while Ca++/Mg++-dependent endonuclease levels increased significantly in relation to day 1. Only at luteolysis (day 17), Ca++/Mg++-dependent endonuclease and MMP-2 spontaneous activity increased significantly. At pregnancy, high levels of MMP-9 and VEGF were observed. These results suggested that during the very early luteal phase, high MMPs activities coupled with high VEGF levels drive the tissue to an angiogenic phenotype, allowing CL growth under LH (Luteinising Hormone) stimulus, while during the late luteal phase, low VEGF and elevate MMPs levels may play a role in the apoptotic tissue and extracellular matrix remodelling during structural luteolysis. In the Experiment III, we described the expression patterns of all distinct VEGF isoforms throughout the oestrous cycle. Furthermore, the mRNA expression and protein levels of both VEGF receptors were also evaluated. Four novel VEGF isoforms (VEGF144, VEGF147, VEGF182, and VEGF164b) were found for the first time in swine and the seven identified isoforms presented four different patterns of expression. All isoforms showed their highest mRNA levels in newly formed CLs (day 1), followed by a decrease during mid-late luteal phase (days 10–17), except for VEGF182, VEGF188 and VEGF144 that showed a differential regulation during late luteal phase (day 14) or at luteolysis (day 17). VEGF protein levels paralleled the most expressed and secreted VEGF120 and VEGF164 isoforms. The VEGF receptors mRNAs showed a different pattern of expression in relation to their ligands, increasing between day 1 and 3 and gradually decreasing during the mid-late luteal phase. The differential regulation of some VEGF isoforms principally during the late luteal phase and luteolysis suggested a specific role of VEGF during tissue remodelling process that occurs either for CL maintenance in case of pregnancy or for noncapillary vessel development essential for tissue removal during structural luteolysis. In summary, our findings allow us to determine relationships among factors involved in the angiogenesis and angioregression mechanisms that take place during the formation and regression of the CL. Thus, CL provides a very interesting model for studying such factors in different fields of the basic research.

Mechanosensitivity in the myenteric plexus of the guinea pig ileum

The enteric nervous system regulates autonomously from the central nervous system all the reflex pathways that control blood flow, motility, water and electrolyte transport and acid secretion. The ability of the gut to function in isolation is one of the most intriguing phenomenons in neurogastroenterology. This requires coding of sensory stimuli by cells in the gut wall. Enteric neurons are prominent candidates to relay mechanosensitivity. Surprisingly, the identity of mechanosensitive neurons in the enteric nervous system as well as the appropriate stimulus modality is unknown despite the evidence that enteric neurons respond to sustained distension. Objectives: The aim of our study was to record from mechanosensitive neurons using physiological stimulus modalities. Identification of sensory neurons is of central importance to understand sensory transmission under normal conditions and in gut diseases associated with sensorimotor dysfunctions, such as Irritable Bowel Syndrome. Only then it will be possible to identify novel targets that help to normalise sensory functions. Methods: We used guinea-pig ileum myenteric plexus preparations and recorded responses of all neurons in a given ganglion with a fast neuroimaging technique based on voltage sensitive dyes. To evoke a mechanical response we used two different kinds of stimuli: firstly we applied a local mechanical distortion of the ganglion surface with von Frey hair. Secondarily we mimic the ganglia deformation during physiological movements of myenteric ganglia in a freely contracting ileal preparation. We were able to reliably and reproducibly mimic this distortion by intraganglionic injections of small volumes of oxygenated and buffered Krebs solution using stimulus parameters that correspond to single contractions. We also performed in every ganglion tested, electrical stimulations to evoke fast excitatory postsynaptic potentials. Immunohistochemistry reactions were done with antibodies against Calbindin and NeuN, considered markers for sensory neurons. Results: Recordings were performed in 46 ganglia from 31 guinea pigs. In every ganglion tested we found from 1 to 21 (from 3% to 62%) responding cells with a median value of 7 (24% of the total number of neurons). The response consisted of an almost instantaneous spike discharge that showed adaptation. The median value of the action potential frequency in the responding neurons was 2.0 Hz, with a recording time of 1255 ms. The spike discharge lasted for 302 ± 231 ms and occurred only during the initial deformation phase. During sustained deformation no spike discharge was observed. The response was reproducible and was a direct activation of the enteric neurons since it remained after synaptic blockade with hexamethonium or ω-conotoxin and after long time perfusion with capsaicin. Muscle tone appears not to be required for activation of mechanosensory neurons. Mechanosensory neurons showed a response to mechanical stimulation related to the stimulus strength. All mechanosensory neurons received fast synaptic inputs. There was no correlation between mechanosensitivity and Calbindin-IR and NeuN-IR (44% of mechanosensitive neurones Calb-IR-/NeuN-IR-). Conclusions: We identified mechanosensitive neurons in the myenteric plexus of the guinea pig ileum which responded to brief deformation. These cells appear to be rapidly accommodating neurons which respond to dynamic change. All mechanosensitive neurons received fast synaptic input suggesting that their activity can be highly modulated by other neurons and hence there is a low stimulus fidelity which allows adjusting the gain in a sensory network. Mechanosensitivity appears to be a common feature of many enteric neurons belonging to different functional classes. This supports the existence of multifunctional enteric neurons which may fulfil sensory, integrative and motor functions.

Outcome respiratorio a lungo termine nei soggetti affetti da cardiopatia congenita sottoposti ad intervento cardiochirurgico

Introduction: In the last years cardiac surgery for congenital heart disease (CHD) reduced dramatically mortality modifying prognosis, but, at the same time, increased morbidity in this patient population. Respiratory and cardiovascular systems are strictly anatomically and functionally connected, so that alterations of pulmonary hemodynamic conditions modify respiratory function. While very short-term alterations of respiratory mechanics after surgery were investigated by many authors, not as much works focused on long-term changes. In these subjects rest respiratory function may be limited by several factor: CHD itself (fetal pulmonary perfusion influences vascular and alveolar development), extracorporeal circulation (CEC), thoracotomy and/or sternotomy, rib and sternal contusions, pleural adhesions and pleural fibrosis, secondary to surgical injury. Moreover inflammatory cascade, triggered by CEC, can cause endothelial damage and compromise gas exchange. Aims: The project was conceived to 1) determine severity of respiratory functional impairement in different CHD undergone to surgical correction/palliation; 2) identify the most and the least CHD involved by pulmonary impairement; 3) find a correlation between a specific hemodynamic condition and functional anomaly, and 4) between rest respiratory function and cardiopulmonary exercise test. Materials and methods: We studied 113 subjects with CHD undergone to surgery, and distinguished by group in accord to pulmonary blood flow (group 0: 28 pts with normal pulmonary flow; group 1: 22 pts with increased flow; group 2: 43 pts with decreased flow; group 3: 20 pts with total cavo-pulmonary anastomosis-TCPC) followed by the Pediatric Cardiology and Cardiac Surgery Unit, and we compare them to 37 age- and sex-matched healthy subjects. In Pediatric Pulmonology Unit all pts performed respiratory function tests (static and dynamic volumes, flow/volume curve, airway resistances-raw- and conductance-gaw-, lung diffusion of CO-DLCO- and DLCO/alveolar volume), and CHD pts the same day had cardiopulmonary test. They all were examined and had allergological tests, and respiratory medical history. Results: restrictive pattern (measured on total lung capacity-TLC- and vital capacity-VC) was in all CHD groups, and up to 45% in group 2 and 3. Comparing all groups, we found a significant difference in TLC between healthy and group 2 (p=0.001) and 3 (p=0.004), and in VC between group 2 and healthy (p=0.001) and group 1(p=0.034). Inspiratory capacity (IC) was decreased in group 2 related to healthy (p<0.001) and group 1 (p=0.037). We showed a direct correlation between TLC and VC with age at surgery (p=0.01) and inverse with number of surgical interventions (p=0.03). Reduced FEV1/FVC ratio, Gaw and increased Raw were mostly present in group 3. DLCO was impaired in all groups, but up to 80% in group 3 and 50% in group 2; when corrected for alveolar volume (DLCO/VA) reduction persisted in group 3 (20%), 2 (6.2%) and 0 (7.1%). Exercise test was impaired in all groups: VO2max and VE markedly reduced in all but especially in group 3, and VE/VCO2 slope, marker of ventilatory response to exercise, is increased (<36) in 62.5% of group 3, where other pts had anyway value>32. Comparing group 3 and 2, the most involved categories, we found difference in VO2max and VE/VCO2 slope (respectively p=0.02 and p<0.0001). We evidenced correlation between rest and exercise tests, especially in group 0 (between VO2max and FVC, FEV1, VC, IC; inverse relation between VE/VCO2slope and FVC, FEV1 and VC), but also in group 1 (VO2max and IC), group 2 (VO2max and FVC and FEV1); never in group 3. Discussion: According with literature, we found a frequent impairment of rest pulmonary function in all groups, but especially in group 2 and 3. Restrictive pattern was the most frequent alteration probably due to compromised pulmonary (vascular and alveolar) development secondary to hypoperfusion in fetal and pre-surgery (and pre-TCPC)life. Parenchymal fibrosis, pleural adhesions and thoracic deformities can add further limitation, as showed by the correlation between group 3 and number of surgical intervention. Exercise tests were limited, particularly in group 3 (complex anatomy and lost of chronotropic response), and we found correlations between rest and exercise tests in all but group 3. We speculate that in this patients hemodynamic exceeds respiratory contribution, though markedly decreased.

Fluidodinamica biomedica computazionale del flusso sanguigno in presenza di affezioni strumentali e/o patologiche

La tesi di Dottorato studia il flusso sanguigno tramite un codice agli elementi finiti (COMSOL Multiphysics). Nell’arteria è presente un catetere Doppler (in posizione concentrica o decentrata rispetto all’asse di simmetria) o di stenosi di varia forma ed estensione. Le arterie sono solidi cilindrici rigidi, elastici o iperelastici. Le arterie hanno diametri di 6 mm, 5 mm, 4 mm e 2 mm. Il flusso ematico è in regime laminare stazionario e transitorio, ed il sangue è un fluido non-Newtoniano di Casson, modificato secondo la formulazione di Gonzales & Moraga. Le analisi numeriche sono realizzate in domini tridimensionali e bidimensionali, in quest’ultimo caso analizzando l’interazione fluido-strutturale. Nei casi tridimensionali, le arterie (simulazioni fluidodinamiche) sono infinitamente rigide: ricavato il campo di pressione si procede quindi all’analisi strutturale, per determinare le variazioni di sezione e la permanenza del disturbo sul flusso. La portata sanguigna è determinata nei casi tridimensionali con catetere individuando tre valori (massimo, minimo e medio); mentre per i casi 2D e tridimensionali con arterie stenotiche la legge di pressione riproduce l’impulso ematico. La mesh è triangolare (2D) o tetraedrica (3D), infittita alla parete ed a valle dell’ostacolo, per catturare le ricircolazioni. Alla tesi sono allegate due appendici, che studiano con codici CFD la trasmissione del calore in microcanali e l’ evaporazione di gocce d’acqua in sistemi non confinati. La fluidodinamica nei microcanali è analoga all’emodinamica nei capillari. Il metodo Euleriano-Lagrangiano (simulazioni dell’evaporazione) schematizza la natura mista del sangue. La parte inerente ai microcanali analizza il transitorio a seguito dell’applicazione di un flusso termico variabile nel tempo, variando velocità in ingresso e dimensioni del microcanale. L’indagine sull’evaporazione di gocce è un’analisi parametrica in 3D, che esamina il peso del singolo parametro (temperatura esterna, diametro iniziale, umidità relativa, velocità iniziale, coefficiente di diffusione) per individuare quello che influenza maggiormente il fenomeno.

Lesion detection and classification in breast cancer: evaluation of approaches based on morphological features, tracer kinetic modelling and semi-quantitative parameters in MR functional imaging (DCE-MRI)

The diagnosis, grading and classification of tumours has benefited considerably from the development of DCE-MRI which is now essential to the adequate clinical management of many tumour types due to its capability in detecting active angiogenesis. Several strategies have been proposed for DCE-MRI evaluation. Visual inspection of contrast agent concentration curves vs time is a very simple yet operator dependent procedure, therefore more objective approaches have been developed in order to facilitate comparison between studies. In so called model free approaches, descriptive or heuristic information extracted from time series raw data have been used for tissue classification. The main issue concerning these schemes is that they have not a direct interpretation in terms of physiological properties of the tissues. On the other hand, model based investigations typically involve compartmental tracer kinetic modelling and pixel-by-pixel estimation of kinetic parameters via non-linear regression applied on region of interests opportunely selected by the physician. This approach has the advantage to provide parameters directly related to the pathophysiological properties of the tissue such as vessel permeability, local regional blood flow, extraction fraction, concentration gradient between plasma and extravascular-extracellular space. Anyway, nonlinear modelling is computational demanding and the accuracy of the estimates can be affected by the signal-to-noise ratio and by the initial solutions. The principal aim of this thesis is investigate the use of semi-quantitative and quantitative parameters for segmentation and classification of breast lesion. The objectives can be subdivided as follow: describe the principal techniques to evaluate time intensity curve in DCE-MRI with focus on kinetic model proposed in literature; to evaluate the influence in parametrization choice for a classic bi-compartmental kinetic models; to evaluate the performance of a method for simultaneous tracer kinetic modelling and pixel classification; to evaluate performance of machine learning techniques training for segmentation and classification of breast lesion.

CT perfusion image processing: analysis of liver tumors

Perfusion CT imaging of the liver has potential to improve evaluation of tumour angiogenesis. Quantitative parameters can be obtained applying mathematical models to Time Attenuation Curve (TAC). However, there are still some difficulties for an accurate quantification of perfusion parameters due, for example, to algorithms employed, to mathematical model, to patient’s weight and cardiac output and to the acquisition system. In this thesis, new parameters and alternative methodologies about liver perfusion CT are presented in order to investigate the cause of variability of this technique. Firstly analysis were made to assess the variability related to the mathematical model used to compute arterial Blood Flow (BFa) values. Results were obtained implementing algorithms based on “ maximum slope method” and “Dual input one compartment model” . Statistical analysis on simulated data demonstrated that the two methods are not interchangeable. Anyway slope method is always applicable in clinical context. Then variability related to TAC processing in the application of slope method is analyzed. Results compared with manual selection allow to identify the best automatic algorithm to compute BFa. The consistency of a Standardized Perfusion Index (SPV) was evaluated and a simplified calibration procedure was proposed. At the end the quantitative value of perfusion map was analyzed. ROI approach and map approach provide related values of BFa and this means that pixel by pixel algorithm give reliable quantitative results. Also in pixel by pixel approach slope method give better results. In conclusion the development of new automatic algorithms for a consistent computation of BFa and the analysis and definition of simplified technique to compute SPV parameter, represent an improvement in the field of liver perfusion CT analysis.

La neoangiogenesi nel reimpianto autologo di tessuto ovarico in pazienti a rischio di fallimento ovarico precoce: ottimizzazione dei metodi di crioconservazione di corticle ovarica

Introduzione L’efficacia dei chemio/radioterapici ha aumentato notevolmente l’aspettativa di vita delle pazienti oncologiche, tuttavia, questi trattamenti possono compromettere la funzionalità ovarica. La crioconservazione di tessuto ovarico, con il successivo reimpianto, ha lo scopo di preservare la fertilità delle pazienti a rischio di fallimento ovarico precoce. Scopo dello studio Definire la migliore procedura di crioconservazione e reimpianto in grado di ottenere la neovascolarizzazione del tessuto reimpiantato nel minor tempo possibile al fine di diminuire la perdita follicolare causata dall’ischemia durante la procedura. Materiali e metodi Per ciascuna paziente (3) le biopsie ovariche, sono state prelevate laparoscopicamente e crioconservate secondo il protocollo di congelamento lento/scongelamento rapido. Campioni di corticale ovarica sono stati processati per l’analisi istologica, ultrastrutturale, immuistochimica e confocale per valutare la preservazione morfologiaca del tessuto. Le fettine di corticale ovarica sono state scongelate e reimpiantate ortotopicamente (2), nelle ovaia e in due tasche peritoneali, o eterotopicamente (1), in due tasche create nel sottocute sovrapubico. Risultati Le analisi di microscopia hanno mostrato il mantenimento di una discreta morfologia dello stroma, e dei vasi criopreservati e un lieve ma non significativo danneggiamento dei follicoli scongelati. Tutte le pazienti hanno mostrato la ripresa della funzionalità endocrina rispettivamente dopo 2/4 mesi dal reimpianto. Il color-doppler, inoltre ha rivelato un significativo aumento della vascolarizzazione ovarica rispetto alla quasi totale assenza di vascolarizzazione prima del reimpianto, quando le pazienti mostravano una conclamata menopausa. Conclusioni Lo studio ha confermato la ripresa della vascolarizzazione dell’ovaio in seguito a reimpianto avascolare di fettine di corticale, senza l’impiego di fattori esogeni o meccanici aggiuntivi, in tempi concordanti con i dati della letteratura. I risultati sono incoraggianti e l’avanzare degli studi e della ricerca potranno contribuire allo sviluppo di nuove metodologie di reimpianto che possano avere un successo clinico ed una sicurezza superiori a quelle finora ottenute.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Aspetti Ultrasonografici in Corso di Sviluppo Placentare e Adattamento Cardiovascolare nella Cavalla Gravida: ecografia bidimensionale, ecocardiografia ed ecocontrastografia (CEUS)

In corso di gravidanza normale avvengono modificazioni emodinamiche centrali e periferiche volte a garantire le crescenti richieste nutritive dell'unità feto-placentare. L’ecografia con mezzo di contrasto (CEUS-Contrast Enhanced Ultrasonography) a base di microbolle offre una nuova opportunità di monitorare e quantificare la perfusione utero-placentare in condizioni normali e patologiche. L’ecocardiografia è stata ampiamente usata in medicina umana per valutare l’adattamento morfo-funzionale cardiaco materno durante la gravidanza. Gli scopi di questo lavoro prospettico sono stati di applicare, per la prima volta nella specie equina, un mezzo di contrasto di II generazione (Sonovue®), al fine quantificare la perfusione utero-placentare in corso di gravidanza normale, valutandone gli effetti sul benessere materno-fetale e di descrivere le modificazioni nei parametri ecocardiografici morfometrici e funzionali cardiaci, in particolare relativi alla funzione del ventricolo sinistro nel corso di una gravidanza fisiologica. Due fattrici sane di razza Trottatore sono state monitorate ecograficamente in maniera seriale durante l’intero corso della gravidanza, tramite esame bidimensionale, ecocontrastografia dell'unità utero-placentare, flussimetria Doppler delle arterie uterine, ecocardiografia materna in modalità bidimensionale, M-mode, Doppler e Tissue Doppler Imaging. I neonati sono stati clinicamente monitorati e gli invogli fetali esaminati. Il pattern di microperfusione utero-placentare è valutabile quali-quantitativamente tramite la CEUS e dimostra un’aumento del flusso a livello di microvascolarizzazione uterina con l'avanzare della gravidanza; non è stata rilevata la presenza di microbolle a livello di strutture fetali nè effetti dannosi sul benessere materno-fetale. In questo studio sono state osservate delle modificazioni cardiache materne in corso di gravidanza fisiologica, relative all'aumento della FC, del CO ed in particolare all'aumento delle dimensioni dell'atrio sinistro ed a modificazioni nelle onde di velocità di flusso e tissutali di riempimento del ventricolo sinistro.

Analisi dei fattori che influenzano la crescita delle arterie polmonari nei pazienti con cuore funzionalmente univentricolare sottoposti ad anastomosi cavo-polmonare bidirezionale

La BCPA e’ ormai utilizzata routinariamente come stadio intermedio della palliazione di Fontan nel trattamento dei cuori funzionalmente univentricolari e si conferma intervento a ridotta mortalità e morbilità. La crescita delle arterie polmonari dopo BCPA è abbastanza variabile ed imprevedibile, le dimensioni indicizzate dei rami principali sembrano tendenzialmente ridotte, tranne quelle dell’arteria lobare inferiore destra, la quale beneficerebbe del flusso preferenziale della BCPA. Proprio i pazienti con arterie polmonari più piccole mostrerebbero un maggior incremento di dimensioni delle stesse dopo BPCA. Il ruolo ed i vantaggi del flusso accessorio nella BCPA sono ancora da definire; tuttavia sembrerebbe offrire vantaggi in termini di crescita delle arterie polmonari, soprattutto il ramo lobare sinistro, e di miglior outcome dopo TCPC.